Static structure factor for graphene in a magnetic field 
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A close study is made of the static structure factor for graphene in a magnetic field at integer 
filling factors v, with focus on revealing possible signatures of "relativistic" quantum field theory 
in the low-energy physics of graphene. It is pointed out, in particular, that for graphene even the 
vacuum state has a nonzero density spectral weight, which, together with the structure factor for all 
u, grows significantly with increasing wave vector; such unusual features of density correlations are 
a "relativistic" effect deriving from massless Dirac quasiparticles in graphene. Remarkably it turns 
out that the zero-energy Landau levels of electrons or holes, characteristic to graphene, remain 
indistinguishable in density response from the vacuum state, although they are distinct in Hall 
conductance. 

PACS numbers: 73.43.-f,71.10.Pm,77.22.Ch 



I. INTRODUCTION 

A great deal of attention has recently been directed 
to graphene, a monolayer of carbon atoms, both ex- 
perimentallyi'^'^ and theoreticallj«^i^iii^i^ii2,. Graphene 
is marked with its novel charge carriers that behave 
like massless Dirac fermions with effective speed of light 
~ 10^ m/s « c/300. It thus provides a special op- 
portunity to study "relativistic" quantum dynamics in 
condensed-matter systems. Experiments have revealed a 
number of exotic transport properties of graphene, such 
as the half-integer quantum Hall (QH) effect and minimal 
conductivity. 

The dynamics of Dirac fermions becomes particularly 
interesting in a magnetic field and leads to peculiar quan- 
tum phenomena, such as fermion number fractionaliza- 
tion^^ and spectral asymmetry ; ^^i^^i^** intimately tied to 
the chiral anomaly in dimensions. Actually the half- 
integer QH effect and the presence of the zero-energy 
Landau levels observedi*^ in graphene are a manifesta- 
tion of fermion number fractionalization. 

It would be important to explore further possible sig- 
natures of relativistic quantum field theory in the low- 
energy physics of graphene. An interesting proposali^ 
along this direction is to simulate the Klein paradox^^ (or 
tunneling) in graphene. Calculationsii of the dielectric 
function also reveal that the electromagnetic response 
of graphene is substantially different from that of con- 
ventional two-dimensional systems. The difference be- 
comes even prominent under a strong magnetic fieldj^^ In 
particular, for graphene the vacuum state is a dielectric 
medium and carries an appreciable amount of electric and 
magnetic susceptibilities over all range of wavelengths; 
this reflects the presence of the "Dirac sea". Curiously 
the zero-energy Landau levels, though distinct in Hall 
conductance, hardly contribute to the susceptibilities. 

The purpose of this paper is to study further aspects 
of the response of graphene in a magnetic field at integer 
filling factors v, with focus on the static structure factor 
s(p) ~ (p-pPp), which is directly related to the cross sec- 
tion in inelastic light scattering by graphene. In partic- 



ular, we point out that graphene has unusual character- 
istics of electronic correlations at short distances, which 
derive from the "relativistic" nature of massless quasipar- 
ticles. It is also shown that the zero-energy Landau levels 
remain indistinguishable in spectral weight (p_pPp) from 
the vacuum state. 

In Sec. II we briefly review the low-energy effective 
theory of graphene in a magnetic field and derive, as 
a preliminary, the structure factor for conventional QH 
systems. In Sec. Ill we study the structure factor for 
graphene. Section IV is devoted to a summary and dis- 



II. LOW-ENERGY EFFECTIVE THEORY 

Graphene has a honeycomb lattice consisting of two 
triangle sublattices of carbon atoms with one electron 
per site. It is a gapless semiconductor and its low-energy 
electronic transport is described by an effective Hamilto- 
nian of the formiS 



H 



n± = VFiaiIli + a2ll2±ma3)-eAa, (2.1) 

where H.; ~ —idi + eAi [i = (1,2) or {x,y)] includes 
coupling to external electromagnetic potentials — 
[Ai^Ao)] vp ~ 10^ m/s is the Fermi velocity. The two- 
component spinors = V'2)* and x = (xii X2)* stand 
for the electron fields near the two inequivalent Fermi 
points {K and K') where the spectrum becomes linear; 
("017X2) reside on the same sublattice and (■02, Xi) on 
another. 

For generality we have introduced a tiny "mass" gap 
TO > which spoils the valley SU(2) symmetry of H. 
Actually the observed 1/ — ±1 Hall plateaus'^ suggest 
such a tiny mass gap4ii2ii£ We keep to ^ to clarify the 
particle- hole character of the lowest Landau levels, but 
practically set to — > 0. 

We suppress the electron spin, which is treated as a 
global SU(2) symmetry of H by doubling the fields and 
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X- The Zeeman splitting, though ignored for simphcity, 
is readily incorporated. 

The Coulomb interaction is written as 



H 



Coul 



p P— p Pp) 



(2.2) 



where pp is the Fourier transform of the electron num- 
ber density p = ip^ip + X^X! ^'p = 27ra/(eb|p|) is the 
Coulomb potential with the fine-structure constant a = 
e^/(47reo) « 1/137 and the substrate dielectric constant 
Ch- We shall discuss the effect of i/Coui ja^g,-^ 

Let us place graphene in a strong magnetic field and 
study how the electrons in graphene respond to a weak 
potential Aq{x). We set Ai B{—y,0) to supply a 
uniform magnetic field = B > normal to the sample 
plane. 

When Aq = 0, the eigenmodes of H are Landau levels 
of tp and X of energy 



2^2/2, 



(2.3) 



labeled by integers n — 0, ±1, ±2, . . ., and Px (or yo = 
with the magnetic length £ = 1/ V eB) ; uic = vy/(. 
is the basic cyclotron frequency. Here s„ = sgn{n} = ±1 
specifies the sign of the energy e„. 

For n ^ 0, V' snd x have the same spectrum sym- 
metric about e = 0. The n — level of i> has negative 
energy eo_ = —v-pm while that of x has positive energy 
£0+ = v-pni; these n = Oip levels represent holes and elec- 
trons via quantization. With the electron spin taken into 
account, each Landau level is thus four-fold degenerate, 
except for the doubly-degenerate n = 0± levels. The 
n = O-t eigenmodes have components only on each sepa- 
rate sublattice. 

To make this Landau- level structure explicit, it is use- 
ful to pass to the |n, yo) basis, with the expansionS£ 

V'(x,t) = En,yo(^l"'2/o)'^"(2/0'*)- (From now on, we 
shall only display the tp sector since the x sector is ob- 
tained by reversing the sign of m.) The Hamiltonian H 
thereby is rewritten as 



H ^ 



n— — oo 



(2.4) 



and the charge density p-p{t) — J (Pxe^'^'^ ^p^■^p as 
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are defined in terms of the harmonic oscillator eigenstates 
{\n)} with a^a|n) — n\n) and [a, a^] — 1; p = Py+ ip^ 
and ~ Pii — iPx ■ More explicitly,— 




kl VV2^ " 



i^2p2] (2.8) 



for k > n, and /„fc(p) = [//cn(-p)]^- Actually /fe„(p) 
are the coefficient functions that characterize the charge 
density for the " nonrelativistic" Hall electrons 



P-P = e 



oo 



^pV4 



k,n=0 



fknip) / dyoV'ie*P-'"^„ (2.9) 



expressed in terms of their eigenmodes {4'n{yoit)} with 
energy e„ = {eB/m*){n + 1/2). 

Let us first consider, as an exercise, the static structure 
factor or the spectral weight {p-p Pp) = {G\p-p Pp\G) for 
a state |G) of free "nonrelativistic" Hall electrons with in- 
teger filling factor v (per spin). One may use Eq. (j2.9|) 
and note, in taking the matrix element {G\p-pPp\G), 

that 'ipk'{y'o,t)ipl{yo,t) dk'kS{yQ - yo) for unoccupied 
levels {k',k) so that the result is proportional to the 
number of electrons per occupied level J dyo tp^^^'n 
LxLy/{2TTi^). (Note here that S{yQ 0) = Lx/{2iTi'^) 
with Lx = J dx.) This yields the spectral weight 

^ oo — 1 

ip-ppp) - A ^^^'^'/^ E E i/^«(p)i' (2-10) 



' 71=0 



for p 7^ 0, where fl — LxLy denotes the total area. 
Let us recall that the static structure factor— 

s{p) = {{p-pPp)-N^Sp^o)/N. (2.11) 

is defined by (p-p Pp) with its p = component isolated, 
where = i/ fl / {2Tr£'^) stands for the total electron num- 
ber. Accordingly, in general, s(p — > 0) = owing to 
charge conservation. Noting the formula 



fe=0 



l//cn(p)|' 



(2.12) 



allows one to cast Eq. (|2.10p into the structure factor at 
filling factor ly 



P-P 



-t 



E gtniP) dyotPle^^-'^t/M, (2.5) .(p) ^ 1 - _ ^ ^ | j^,^(p)|2^ (2.13) 

n — n I — n 



fc,n— — oo 

with the coefficient matrix 

5nn'(P) = 5 C+C+ /|„|_ij„,|_i(p) 



which agrees with a known result. In particular 



n=0 fc=0 
23 



.(p)|..i = l-e-^V/2. 



(2.14) 



+s„s„.c„ c„, /|„|,|„'|(p) (2.6) 



and = ^/1± VFm/en, r = (ri,r2) = {W^d / dyo,Vo) 
stands for the center coordinate with uncertainty 
[^i,''2] — ■ Here the coefficient functions 



Note first that s(p ^ 0) — > since fkn{0) = Skn- Note 
also that s(p) — > 1 as p — > oo for all v. To see what 
this means let us recall the following: For a collection of 
classical particles, s(p) is written as^ 



/fcn(p) 



|g-i(p/V2)at g-i(pt/V2)a 



In) (2.7) 



sip) - 1 = 



Ar 



(2.15) 
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for p ^ 0, i.e., as an average over the relative positions 
Ar of particles surrounding a given particle. As a result, 
s(p) ^ 1 for p ^ oo if Ar ^ 0, e.g., for particles formu- 
lated on a lattice. Such behavior of s(p) thus implies the 
absence of particle correlations at short distances. 

III. STATIC RESPONSE OF GRAPHENE 

In this section we study the case of graphene. In the 
present treatment the t/i and x sectors are independent 
and the spectral weight (pp) with p ~ + given 
by the sum {p^ p^) + {p^p^) for p 7^ 0. For the -0 sector 
one may simply replace, in Eq. (|2.10p . /fc„(p) by gf^iv) 
of Eq. (|2.6[) and note that the level indices (fc, n) now run 
over all integers, 0_, ±1, ±2, • • •. Analogously, (7fc„(p) for 
the X sector is obtained from .g^„(p) by setting m —m 
(i.e., c+ ^ c~), or equivalently 

5L(P)=A-«(P)- (3-1) 

Let us denote V^iz) = |gt(p)P and T^z) = |gL(p)P 
for short. They actually are functions of z = ^^p^/2, are 
thus symmetric in {k,n), and have the property 

Tfcn = r^fc__„, Tfcn = r^fc, (3.2) 
with i = -tjj, X- The spectral weight is now written as 

{p-pPp)/n = {Xs/2TTf) e-i'"p\l^ + F), (3.3) 

where P ~ Yl'k'l2n^kni^)' K — stands for the spin 
degeneracy. One can write as 

k>j n<j 

when the ip Landau levels are occupied up to the jth level 
(j = 0_, ±1, . . .); analogously for 1^-. 

Note that Eq. relates I''' and so that 



The equalities = I^^ and // + If = ^^(^+1) + 
^^(j+i) then imply the following: (i) The ip and x sec- 
tors contribute equally to the vacuum spectral weight, 
{P-pPp)\u=o « Iq- + ^-1 = 2-^o'_- (ii) The spectral 
weight (p-pPp) is the same for the charge-conjugate 
states with i/ — ±2(2j -1-1). In view of this we shall 
focus on the case v > from now on. 

In the limit m — > of practical interest, distinctions 
between F'^ and F^ disappear: F^„ = F^„ = ^-k.-n fo^' 
kn ^ and F^q_ = ^±ko+ ^'^^ 7^ 0, as seen from 
Eq. (EH). This yields 




1 2 3 4 5 6 



(b) 




FIG. 1: (a) Profile of the transition rate e-^Vf^iz) for 
m ^ 0, labeled with {k,n). (b) Vacuum spectral weight 
(p-pPp)jiy=o for N = 800 (dashed line) and the difference 
— 5(P-pPp)|i^=o = {pp)^^^ — (pp)^^" (real line), both in units 
of Kn/{2-Kf). 



In particular, /q_ = /g^ = I^^ = shows that the 
graphene vacuum and the v = ±2 states, all of zero 
energy, have the same spectral weight (p-pPp)|i/=2 cx; 

+ /q^ — 2Iq^ . Actually we earlier noted such degen- 
eracy in electric susceptibility for the 1^ = and i/ = ±2 
stateSfi^ and these states now turn out to be indistin- 
guishable in both the real and imaginary parts of density 
response; we discuss this point in more detail later. 

To proceed further let us note some basic properties 
of the transition rates F^„(z) with z — ^^p^/2. It is 
clear from Eq. (|3.2p that the transition rates among 
the positive-energy states and those among the negative- 
energy states are essentially the same whereas they are 
different from the rates across the Dirac sea (ex Ffe„ with 
kn < 0). This feature is seen from Fig. 1 (a), which il- 
lustrates typical profiles of e~^rfe„(2:) for m ^ 0. In 
general, e~^F^„(z) is significantly peaked at some value 
of i\p\ and the peak position rises only gradually as the 
level gap ~ \k ~ n\ increases. The structure factor s(p) 
for small |p| is thus governed by virtual transitions to 
neighboring levels while its property at larger |p| is de- 
termined by transitions across larger gaps. 

It is readily expected from this property that the sum 
over an infinite number of negative-energy levels would 
make the weight (p_pPp) cutoff (N) dependent for small 
p^, actually to 0{£^p^). In particular, the vacuum weight 

{P-pPp)l=o (X 2e"^Ef=iEiLort-n(^) vanishes at 
z = and grows rapidly with z for fixed TV, whereas 
it diverges as ^ 00 for z 7^ 0. See Fig. 1 (b). In- 
deed, evaluating the 0{z) terms in Tk,-n{z) shows that 
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the divergence is logarithmic in N, 



N N 
fc=ln=0 



N-1 

iz5:(^/fcTT-^/fc)^ 

k=0 

« (l/8)zln(ciiV) 



(3.7) 



with ci « 54.088 and for m ^ 0. This ultraviolet diver- 
gence is physical. The infinite depth of the Dirac sea, of 
course, is an artifact of the continuum model (j2.1[) and 
the cutoff scale lOcV^ is set by the energy scale above 
which the model loses its validity. 

Let us recall here that for conventional QH systems 
the charge operator trivially annihilates the vacuum, 
pliy = 0) = 0, and the vacuum spectral weight vanishes. 
Accordingly, the nonzero vacuum weight {p^pPp)\^^o it- 
self is a " relativistic" signature of graphene. 



graphene 



(3.8) 



which is a consequence of particle-hole pair creation or 
the presence of the Dirac sea. 

It is perfectly legitimate to consider the vacuum weight 
with such a physical cutoff N (apart from its precise 
value). One can equally well extract cutoff-insensitive 
information out of it. One possible way is to consider a 
variation of {p-pPp)\u=o for B ^ and B ~ 0. Exper- 
imentally this means measuring the vacuum weight for 
B and B = separately. 

For B = the spectral weight for the graphene vacuum 
|0) with TO ^ is written as 



1 



\p-pPp\0) 



B=0 



d^k 1 

(2^2 



cos 0k+p,k), 

(3.9) 

which represents a collection of virtual transitions^^ from 
a negative-energy electron state with momentum k to a 
positive-energy state with k -I- p via the charge density 
p; here 0k+p,k denotes the angle between k and k -I- p; 
As = 2 and Ay = 2 denote the spin and valley degeneracy. 
This is again ultraviolet divergent. For regularization we 
cut off the k integral at |k| = A and choose the "Fermi 
momentum" A so that the Dirac sea accommodates the 
same number of electrons as in the B case, iVsca = 
AsAvAV(47r) « AsAv(iV + l/2)/(27r£2), i.e., A^ « 2N/e^. 
A direct calculation yields 



(P-pPp) 



3=0 



p2 16A2 
A,Av:g-log^. (3.10) 
327r p^ 



The cutoff (N) dependence thus correctly disappears 
from the difference S{pp) = (pp)^^^ — with the 

result 

6{p.pPp)\,=o/^ - iXs/27:i^)6I{z), 

SI{z) = 2e-'-lt_{z) - (z/4) log(16A^/z). (3.11) 

As seen from Fig. 1 (b), the vacuum weight {p-pPp)\v=a 
differs only slightly for _B 7^ and _B = 0. 



Actually for all v the spectral weight (p_pPp) contains 
the vacuum fluctuations and the static structure factor 
s(p) is necessarily cutoff-dependent. A possible cutoff- 
independent measure in experiment is to compare the 
spectral weights for and v = (with B ^ 0). 

Correspondingly let us define by 



As(p) = {{p-pPp) - (p-pPp)k=o}/A^o 



(3.12) 



the structure factor with the vacuum contribution sub- 
tracted. For the v = 2{2j + 1) states (with j — 0,1, . . .), 
As(p) e-''{Alf + AIf )/{2j + l) in terms of the cutoff- 



independent deviations A/, 



■4, 



-'o„ 



and Alj 



Ij- — I^i- For j > they are rewritten as 



N j 

E E ^k,n 
k=j+l n=0+ 



j N 

EE-'^ 

fc=0+ n=l 



(3.13) 



E^"( 



z e 



EEr^,»(-); (3-14) 

fc=l n=0+ 



analogously for A/^ . In reaching the second line, we 
have used the formulas (valid for m — > and N ~^ 00): 

^ 1 

Er±M(^) = ^(e^-i), 



fe=i 

N 



fe=i 

where z = i?^p^/2 and 

Fniz) = 



1 zl"! 

4NT'(^-1^) 



JV 



fe=i 



Re[fn-iM-ii-p) fkAp)h (3.16) 



Fn{0) = 1 and Fn{z) ^ as z ^ 00; the decrease is 
slower for larger n, as depicted in Fig. 2 (a). 

The subtracted structure factor is now written as 



As(p) 



2j 



^[E^^^(^) 

n— 1 



k = l Tl = + 

(3.17) 

for ly = 2{2j + 1). In particular, As(p)|,y=2 = and 



As(p)|j,=6 



2 r 
3 



^^i(^) 



1 



- z + - z^ 
2 4 



(3.18) 



In Fig. 2 (b) and (c) we plot As(p) for v — 6, 10, 14 and 
for higher filling — 22,42,62. There are some notable 
features: (1) As(p) for |p| 0. This is a conse- 
quence of the uniformity of charge at long wavelengths. 
(2) As(p) for IpI 00. This shows that the large 
IpI portion of s(p) is common to all v and is governed 
by the vacuum weight {p-pPp)\^=o- As a result, s(p) 
itself would rise with |p|, in sharp contrast to the stan- 
dard behavior s(p) ^ 1 as |p| ^ 00 for " nonrelativistic" 
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2 4 6 8 \p\e 10 




2 4 6 8 |p|£ 10 



FIG. 2; (a) Fn{z). (b) Subtracted static structure factor 
As(p) for u^6, 10, 14. (c) As(p) for i/=22, 42, 62. 




2 4 6 8 |p|£ 10 



FIG. 3: (a) Structure factor s'^'(p restricted to the positive- 
energy sector, for u — 2 ^ 10. (b) As(p) a.t u = 2 develops 
as the mass gap becomes large vpm/u!c = 0.1, 1, 10, 100. 

Hall electrons in Eq. (|2.13|) . (3) As(p) tends to saturate 
around 0.25 over a certain broad range of |p| for higher 
fiUing factors v. The subtracted structure factor As(p), 
unlike s(p), refers to quantum fluctuations of a finite 
number of filled levels above the vacuum and is governed 
by virtual transitions, newly allowed or suppressed by the 
presence of such levels, as seen from Eq. p. 131) . Satura- 



tion is a consequence of a competition between the tran- 
sitions among positive-energy states and the suppressed 
vacuum polarization effect. 

In this connection, let us try to retain only the positive- 
energy states and transitions among them. In Fig. 3 (a) 
we plot a projected structure factor s(+)(p) calculated 
from the first term on the right-hand side of Eq. (|3.13p . 
Note that 

s(+)(p)^l/2 for |p|->oo; (3.19) 

i.e., even the positive-energy sector alone does not re- 
cover the s(p) — > 1 behavior of conventional QH systems. 
While s'-^'(p) itself is not directly observable in exper- 
iment, this feature (|3.19p would indicate indirectly that 
As(p) is significantly smaller than 1/2. 

Some remarks on the influence of a mass gap are in or- 
der here. The relativistic treatment makes sense when 
the mass gap is tiny mi <C 1. In the " nonrelativis- 
tic" limit ml 3> ViV where the mass gap is large com- 
pared with the depth of the Dirac sea, in contrast, the 
Landau-level sum is limited to a finite interval N ^ to^£^ 
and ceases to yield a divergence. The virtual transition 
rates across the mass gap scale like Tk^-n oc l/(rn€)^ 
(while Tkn ^ 0{1) for kn > 0), and the vacuum spec- 
tral weight (p_pPp)|;,=o fx J2k,n^':^-n tends to zero like 
{m£)^'^ \og{m'^P) ^ as m£ — > oo. Figure 3 (b) shows 
the manner how As(p)|y=2 develops as the mass gap gets 
large for v-pm/ujc = 0.1 — > 100. The v = 2 state is now 
distinguishable in s(p) from the vacuum. In the ml —t oo 
limit, As(p) w s(p) w s'-^Hp) recovers the standard be- 
havior s(p) — > 1 as IpI — > oo for all integer i^. 

To explore the origin of the unconventional behav- 
ior (|3.19p let us note the formula 

oo 

E l5;c„(p)P = e^V/2. (3 20) 

A;— — oo 

which follows from Eq. (|3.15p . Actually this formula is 
valid for m 7^ as well, since it is independently derived 
from the evaluation of the static weight {G\p-pPp\G) for 
a hypothetical state \G) with only one filled Landau level 
{n} and all other levels empty. It is thus a consequence 
of completeness of physical states. 

In the " nonrelativistic" limit the Dirac sea effectively 
disappears. Then Eq. ([X^ is reduced to Eq. ((TT^ . 
which in turn leads to the behavior s(p) 1 for |p| 
00. 

For TO — > the large- |p| behavior of Eq. p.20p is gov- 
erned by the sums over \k\ S> |n| so that 

00 

e-''p'/'El-9±fc-(P)l'^ 1/2, (3.21) 

k~0 

as seen also from Eq. (|3.15p . This precisely accounts for 
the behavior of s(+'(p) in Eq. (|3.19p . It is now clear that 
both the masslessness (or a tiny mass gap) of the quasi- 
particles and the presence of the Dirac sea are crucial for 
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the unusual density correlations s'^+^(p) 1/2 at short 
distances. In other words, correlations fail to vanish at 
short distances because of particle-hole pair creation; the 
inability of localizing massless particles is a purely "rel- 
ativistic" effect, underlying also the Klein paradoxji^ 

Equation (I3.20|) may suggest that the vacuum spec- 
tral weight (p-pPp)|iy=o would have a contribution of 
0{N) from the Dirac sea, It actually rises like InA^, as 
we have seen in Eq. (|3.7p . Had we defined the " vacuum" 
static structure factor s(p) with normalization by the to- 
tal number of electrons in the Dirac sea, it would vanish 
like s(p) c>c (lnA^)/iV 0. This shows again that, for 
massless particles and in the presence of the Dirac sea, 
the standard (classical) picture of particle correlations 
does not necessarily make sense. 

Finally we wish to comment on some peculiarities of 
the zero-energy n — 0± levels, especially on how the 
Coulomb interaction Vp — 27rck/(eb|p|) affects their re- 
sponse. In the random-phase approximation (RPA) the 
effect of ff*^""' of Eq. ([^ is readily included^^ into the 
polarization function ('^ —i{pp)), 

Prpa(p, l^) = P(p, u;)/{l - Vp P(p, uj)}, (3.22) 

once one knows the polarization function P(p, lu) for non- 
interacting Hall electrons. Suppose now that we start 
with the vacuum state and fill up the n — 0+ level to 
reach the = 2 state. Then P{p, oj) will change by an 
amount proportional to 



AP(p,c^) cx J2 



fc>i 



E 



rl,-Jp), (3.23) 



fc>i "+ *- 



which is easily seen to vanish for zero mass-gap m — > 0. 
This shows that the v = Q and v — ±2 states remain 
indistinguishable in density response even at the RPA 
level, indicating the robustness of the zero modes against 
perturbations. In particular, the static structure fac- 
tor s(p) ~ / dcj Prpa(p, 1.^) and the dielectric function 
e(p, uj) — 1 — Vp P(p, ijj) remain the same for = and 
V = ±2. 



with full of virtual particle-hole pairs. In other words, 
p\i' = Q) ^ Q for graphene while p = 0) = for stan- 
dard QH systems. Correspondingly the graphene vacuum 
state has a nonzero density spectral weight {p-pPp)\u=f)i 
which actually diverges with the (cutoff) depth of the 
Dirac sea and generally rises with |p| significantly. Ex- 
perimentally the nonzero vacuum weight (p_pPp)|i/=o it- 
self as well as its rise with |p| , measured via inelastic light 
scattering, would be a clear signal of the quantum nature 
of the graphene vacuum state. 

For graphene the static structure factor s(p) grows 
with IpI for all in sharp contrast to the behavior 
s(p) ^ 1 as IpI — > oo of conventional QH systems; this 
is because of pair creation at short distances. Actually, 
even the transitions among positive-energy states fail to 
recover the standard behavior and lead to s^+^(p) ^1/2 
for IpI — + oo, as noted in Eq. (|3.19p . This feature 
of s'^'(p) itself, unfortunately, is not directly observ- 
able. Nevertheless, for large filling factor v, it competes 
with the vacuum polarization effect and would make the 
observable subtracted structure factor As(p) (which is 
sensitive to quantum fluctuations of cutoff-independent 
longer wavelengths) saturate around 1/4 over a certain 
broad range of |p|; this could be a possible signature of 
the underlying dynamics. Such unusual features of par- 
ticle correlations at short distances are a " relativistic" 
effect coming from massless Dirac particles in graphene. 

Of special interest, in addition, are some peculiar fea- 
tures of the lowest {n — Q±) Landau levels of zero energy. 
Remarkably they hardly affect the density response. The 
virtual transitions, both allowed and suppressed anew by 
the presence (or absence) of the n = 0± levels, combine 
to leave the density response unchanged, and this com- 
pensation persists even when the Coulomb interaction is 
taken into account in the RPA. The graphene vacuum 
[v = 0) state and the i' = ±2 states would thus be indis- 
tinguishable in spectral weight (p-pPp) ~ s(p) as well 
as in electric susceptibility, although they are distinct 
in Hall conductance. Experimentally one would expect a 
definite nonzero signal for (p_pPp) which largely remains 
the same over the range |z^| < 2, except for some values of 
v where nontrivial dynamics such as the fractional quan- 
tum Hall effect may come into play. 
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